ASIAC 


REPORT  NO.  1182.10 


ANALYSIS  OF  DAMPED  PLATES 


November  1982 


Approved  for  public  release; 
distribution  unlimited 


AEROSPACE  STRUCTURES 
INFORMATION  AND  ANALYSIS  CENTER 

OPERATED  FOR  THE  AIRFORCE  FUGHT  DYNAMICS  LABORATORY 
BY  ANAMET  LABORATORIES,  INC. 


19990518  068 


Aerospace  Structures 
Information  and  Analysis  Center 


Report  No.  1182. 1C 


This  report  presents  three  design-oriented  methods  for  the 
dynamic  analysis  of  sandwich  plates,  i.e.,  laminated  plates 
incorporating  a  core  layer  of  viscoelastic  material  for  vibration 
damping.  The  methods  are  complementary  in  that  each  represents 
a  different  trade-off  between  generality,  accuracy,  and  cost  of 
use. 

This  work  was  performed  by  the  Aerospace  Structures  Infor¬ 
mation  and  Analysis  Center,  which  is  operated  for  the  Flight 
Dynamics  Laboratory  by  Anamet  Laboratories,  Inc.,  under 
Contract  F33615-81-C-3201 .  The  work  was  performed  in  support 
of  the  Advanced  Radiation  Technology  (ART)  Program  Office 
at  the  Air  Force  Weapons  Laboratory. 

Prepared  by : 


s 

David  A. 

Kienholz , 

(jPh.D. 

Senior  Engineer 


Conor  D.  Jq5?i!hson,  Ph.D. 
Principal  Engineer 


Jatin  Parekh 
Engineer 


Approved  by: 


E.  B.  Flora 
Program  Manager 


Aerospace  Structures 
Information  and  Analysis  Center 

TABLE  OF  CONTENTS 

Page 

1.0  INTRODUCTION  . . .  1 

2 . 0  MODAL  STRAIN  ENERGY  METHOD  .  3 

2.1  OVERVIEW . 3 

2.2  THEORY . 5 

2.3  FINITE  ELEMENT  MODELING  OF  THREE-LAYER  PLATES  .  8 

2.3.1  Choice  of  Elements .  8 

2.3.2  Reduction  of  Equations  of  Motion  ....  10 

2.3.3  Solution  Method . 11 

2.4  EXAMPLE . 13 

3.0  DESIGN  CHARTS  FOR  SANDWICH  PLATES  .  23 

3 . 1  DATA  FORMAT . 2  3 

3.2  DESIGN  CHARTS . 27 

3.2.1  PTU  Boundary  Conditions  . . 27 

3.2.2  PLR  Boundary  Conditions . 61 

3.2.3  PTR  Boundary  Conditions . 71 

3.2.4  PWU  Boundary  Conditions . 81 

3.2.5  PWR  Boundary  Conditions . 92 

3.3  EXAMPLE . 102 

4.0  CLOSED  FORM  SOLUTION  FOR  HIGH-ORDER  MODES  .  106 

4.1  THEORY . 106 

4.2  SOFTWARE  IMPLEMENTATION . .  •  109 

4.3  EXAMPLE . HO 

5.0  SUMMARY  AND  CONCLUSIONS . 12  2 

REFERENCES . 12  4 

APPENDIX  A  -  SAMPLE  INPUT  AND  OUTPUT  FOR 

NASTRAN/MODAL  STRAIN  ENERGY  ANALYSIS 

OF  A  SANDWICH  PLATE  . . A-1 


Aerospace  Structures 
Information  and  Analysis  Center 


LIST.  OF.  FIGURES 


Figure  No. 


Page  No . 


1 

2 


3 


4 


5 


6 

7 

8 

9 


10 


Finite  element  modeling  of  a  sandwich 

panel  with  viscoelastic  core  . . 

Damping  of  the  first  mode  of  simply 

supported  rectangular  sandwich  plate 

obtained  by  NASTRAN/Modal  Strain  Energy 

method  and  by  exact  complex  eigenvalue 

solution  [8]  . 

Damping  of  the  second  mode  of  simply 

supported  rectangular  sandwich  plate 

obtained  by  NASTRAN/Modal  Strain  Energy 

method  and  by  exact  complex  eigenvalue 

solution  [8]  . 

Damping  of  the  third  mode  of  simply 

supported  rectangular  sandwich  plate 

obtained  by  NASTRAN/Modal  Strain  Energy 

method  and  by  exact  complex  eigenvalue 

solution  [8]  . 

Damping  of  the  fourth  mode  of  simply 

supported  rectangular  sandwich  plate 

obtained  by  NASTRAN/Modal  Strain  Energy 

method  and  by  exact  complex  eigenvalue 

solution  [8] . .  . . . 

Damping  of  a  sandwich  rectangular  plate, 

PTU  boundary  conditions , Axy  =1.1,Y=3.5  .  .  .  . 

Damping  of  a  sandwich  rectangular  plate, 

PTU  boundary  conditions , Axy  =2.0,Y=3.5  .  .  .  . 

Damping  of  a  sandwich  rectangular  plate, 

PTU  boundary  conditions,  Axy  =4.0,Y=3.5  .  .  . 

Damping  of  a  sandwich  rectangular  plate, 

PTU  boundary  conditions , Axy  =  1.1,  modes  1 

and  3 ,  variable  Y  and  g  . 

Damping  of  a  sandwich  rectangular  plate, 

PTU  boundary  conditions , Axy  =  1.1,  modes  2 

and  4 ,  variable  Y  and  g . .  .  . 


9 


15 


16 


17 


18 

28 

29 

30 


31 


32 


Aerospace  Structures 
Information  and  Analysis  Center 


LIST  OF  FIGURES  (Continued) 

Figure  No. 

11  Damping  of  a  sandwich  rectangular  plate, 

PTU  boundary  conditions , Axy  =  2,0,  modes  1 
and  3,  variable  Y  and  g  . 

12  Damping  of  a  sandwich  rectangular  plate, 

PTU  boundary  conditions , Axy  =  2.0,  modes  2 
and  4,  variable  Y  and  g . 

13  Damping  of  a  sandwich  rectangular  plate, 

PTU  boundary  conditions , Axy  =  4.0,  modes  1 
and  3 ,  variable  Y  and  g  . 

14  Damping  of  a  sandwich  rectangular  plate, 

PTU  boundary  conditions , Axy  =  4.0,  modes  2 
and  4 ,  variable  Y  and  g  . 

15  Natural  frequencies  of  sandwich  rectangular 
plate,  PTU  boundary  conditions,  Axy  =  1.1, 

Y  =  0.5 . .  . 

16  Natural  frequencies  of  sandwich  rectangular 
plate,  PTU  boundary  conditions,  Axy  =  2.0, 

Y  =  0.5 . . 

17  Natural  frequencies  of  sandwich  rectangular 
plate,  PTU  boundary  conditions,  Axy  =  4.0, 

Y  =  0.5  . 

18  Natural  frequencies  of  sandwich  rectangular 
plate,  PTU  boundary  conditions,  Axy  =  1.1, 

Y  =  1.5 . 

19  Natural  frequencies  of  sandwich  rectangular 
plate,  PTU  boundary  conditions,  Axy  =  2.0, 

Y  =  1.5 . . . 

20  Natural  frequencies  of  sandwich  rectangular 
plate ,  PTU  boundary  conditions ,  Axy  =  4.0, 

Y  =  1.5  .  .  .  . . . 

21  Natural  frequencies  of  sandwich  plate,  PTU 

boundary  conditions,  Axy  =  1*1»  Y  =  3.5  .  . 

22  Natural  frequencies  of  sandwich  rectangular 

plate,  PTU  boundary  conditions,  Axy  =  2.0, 
Y=2.3  . 


Page  No. 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 


Aerospace  Structures 
Information  and  Analysis  Center 


23  Natural  frequencies  of  sandwich  rectangular 
plate,  PTU  boundary  conditions,  Axy  =  4.0, 

Y  =  3.5  .....  . . 

24  Natural  frequencies  of  sandwich  rectangular 
plate,  PTU  boundary  conditions,  Axy  =  1.1, 

Y  =  4.5  ......  . . 

25  Natural  frequencies  of  sandwich  rectangular 
plate,  PTU  boundary  conditions,  Axy  =  2.o, 

Y  =  4.5  . 

26  Natural  frequencies  of  sandwich  rectangular 
plate,  PTU  boundary  conditions,  Axy  =  4.0, 

Y  =  4.5 . . 

27  Damping  of  a  sandwich  rectangular  plate, 

PLR  boundary  conditions,  Axy  =1.1,Y=3.5  .  .  . 

28  Damping  of  a  sandwich  rectangular  plate, 

PLR  boundary  conditions,  Axy  =2.0,Y=3.5  .  .  . 

29  Damping  of  a  sandwich  rectangular  plate, 

PLR  boundary  conditions,  Axy  =4.0,  Y=  3.5  .  .  . 

30  Natural  frequencies  of  a  rectangular 
sandwich  plate,  PLR  boundary  conditions, 

Axy  =1.1,Y=3.5  . 

31  Natural  frequencies  of  a  rectangular 
sandwich  plate,  PLR  boundary  conditions, 

Axy  =2.0,Y=3.5  .  . 

32  Natural  frequencies  of  a  rectangular 
sandwich  plate,  PLR  boundary  conditions, 

Axy  =4.0,Y=3.5  .  . 

33  Damping  of  a  sandwich  rectangular  plate, 

PTR  boundary  conditions,  Axy  =1.1,Y=3.5  .  .  . 

34  Damping  of  a  sandwich  rectangular  plate, 

PTR  boundary  conditions ,  Axy  =  2 . 0 ,  Y  =  3 . 5  .  .  . 

35  Damping  of  a  sandwich  rectangular  plate, 

PTR  boundary  conditions,  Axy  =4.0,Y=3.5  .  .  . 


Page  No. 

45 

46 

47 

48 

62 

63 

64 

65 

66 

67 

72 

73 

74 


Aerospace  Structures 
Information  and  Analysis  Center 


LIST  OF  FIGURES  (Continued) 


Figure  No. 

36  Natural  frequencies  of  a  rectangular 
sandwich  plate,  PTR  boundary  conditions, 

.  Axy  =  1.1,Y=3.5 . .  •  • 

37  Natural  frequencies  of  a  rectangular 
sandwich  plate,  PTR  boundary  conditions, 

Axy  =2.0,Y=3.5  . 

38  Natural  frequencies  of  a  rectangular 
sandwich  plate,  PTR  boundary  conditions, 

Axy  =4.0,Y=3.5  . 

39  Sandwich  plates  with  PWU  and  PWR 

boundary  conditions  . 

40  Damping  of  a  sandwich  rectangular  plate, 

PWU  boundary  conditions,  Axy  =1.1,Y=3.5  .  . 

41  Damping  of  a  sandwich  rectangular  plate, 

PWU  boundary  conditions,  Axy  =2.0,Y=3.5  .  . 

42  Damping  of  a  sandwich  rectangular  plate, 

PWU  boundary  conditions ,  Axy  =4.0,Y=3.5  .  . 

43  Natural  frequencies  of  a  rectangular 
sandwich  plate,  PWU  boundary  conditions, 

Axy  =1.1,Y=3.5  . 

44  Natural  frequencies  of  a  rectangular 
sandwich  plate,  PWU  boundary  conditions, 

Axy  =2.0,Y=3.5  .  . 

45  Natural  frequencies  of  a  rectangular 
sandwich  plate,  PWU  boundary  conditions, 

Axy  =  4.0,  Y  =  3.5  . 

46  Damping  of  a  sandwich  rectangular  plate, 

PWR  boundary  conditions,  Axy  =1.1,  Y=  3.5 

47  Damping  of  a  sandwich  rectangular  plate, 

PWR  boundary  conditions,  Axy  =2.0,Y=3.5  .  • 

48  Damping  of  a  sandwich  rectangular  plate, 

PWR  boundary  conditions,  Axy  =4.0,Y=3.5  .  . 


Aerospace  Structures 
Information  and  Analysis  Center 


Figure 

49 

50 

■  51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


LIST  OF  FIGURES  (Continued) 


Natural  frequencies  of  a  sandwich 
rectangular  plate,  PWR  boundary 
conditions,  Axy  -  1.1,  Y  =  3.5  . 

Natural  frequencies  of  a  sandwich 
rectangular  plate,  PWR  boundary 
conditions,  Axy  =2.0,  Y  =  3.5  .  • 

Natural  frequencies  of  a  sandwich 
rectangular  plate ,  PWR  boundary 
conditions,  Axy  =  4.0,  Y  =  3.5  . 


Effect  of  boundary  conditions  on  damping 
for  rectangular  sandwich  plates , 
g  =  1,3  and  10  . 

Effect  of  boundary  conditions  on  damping 
for  rectangular  sandwich  plates,  g  =  40  . 


Effect  of  boundary  conditions  on  damping 
for  rectangular  sandwich  plates,  g  =  100  . 

Effect  of  boundary  conditions  on  damping 
for  sandwich  plates,  g  =  300  and  1000  .  . 

Effect  of  boundary  conditions  on  damping 
for  rectangular  sandwich  plates, 
g  =  1,3  and  10  .  .  . . 

Effect  of  boundary  conditions  on  damping 
for  rectangular  sandwich  plates,  g  =  40 

Effect  of  boundary  conditions  on  damping 
for  rectangular  sandwich  plates,  g  =  100  . 


Effect  of  boundary  conditions  on  damping 
for  rectangular  sandwich  plates,  g  =  300 

Effect  of  boundary  conditions  on  damping 
for  rectangular  sandwich  plates,  g  =  1000 


Page  No . 


96 

97 

98 

112 

113 

114 

115 

116 

117 

118 

119 

120 


Aerospace  Structures 
Information  and  Analysis  Center 


LIST  OF  TABLES 


Table  No. 


1 


2 


3 


4 


5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


Modal  Frequencies  and  Modal  Loss  Factors 
for  a  Rectangular  Sandwich  Plate, 

Axy  =  1.1,  Y  =  3.5,  =  0.1  . 

Modal  Frequencies  and  Modal  Loss  Factors 
for  a  Rectangular  Sandwich  Plate , 

Axy  =  1.1,  Y  =  3. 5,  =  0. 3  . 

Modal  Frequencies  and  Modal  Loss  Factors 
for  a  Rectangular  Sandwich  Plate, 

Axy  =  1.1,  Y  =  3.5,  q^  =  1.0  . 

Modal  Frequencies  and  Modal  Loss  Factors 
for  a  Rectangular  Sandwich  Plate, 

Axy  =1.1,Y=3.5,  q^=2.0  . 

Reference  Frequencies  . 

Modal  Frequencies  and  Modal  Loss  Factors , 
PTU,  Axy  =1.1,Y=0.5  . 

Modal  Frequencies  and  Modal  Loss  Factors, 
PTU,  Axy  =1.1,Y=1.5  . 

Modal  Frequencies  and  Modal  Loss  Factors, 
PTU,  Axy  =l.l,Y-3.5  ...  . 

Modal  Frequencies  and  Modal  Loss  Factors , 
PTU,  Axy  =1.1,Y=4.5  . 

Modal  Frequencies  and  Modal  Loss  Factors, 
PTU,  Axy  =2.0,Y=0.5  ........ 

Modal  Frequencies  and  Modal  Loss  Factors, 
PTU,  Axy  =2.0,Y=1.5  . 

Modal  Frequencies  and  Modal  Loss  Factors , 
PTU,  Axy  =2.0,Y=3.5  . 

Modal  Frequencies  and  Modal  Loss  Factors, 
PTU,  Axy  =2.0,Y=4.5  . 

Modal  Frequencies  and  Modal  Loss  Factors, 
PTU,  Axy  =4.0,Y=0.5  . 

Modal  Frequencies  and  Modal  Loss  Factors, 
PTU,  Axy  =  4.0,  Y  =  1.5  .  .  .  .  .  .  .  . 


Page  No. 


19 


20 


21 


22 

26 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 


Aerospace  Structures 
Information  and  Analysis  Center 


LIST  OF  TABLES  (Continued) 

Table  No. 


16 

Modal  Frequencies 
PTU,  Axy  =  4.0,  Y 

and 
=  3. 

Modal 

5  .  . 

Loss 

•  •  • 

Factors , 

♦  •  •  •  ® 

o  e 

• 

59 

17 

Modal  Frequencies 
PTU,  Axy  =  4.0,  Y 

and 
=  4. 

Modal 

5  .  . 

Loss 

Factors , 

•  •  •  •  • 

o  • 

• 

60 

18 

Modal  Frequencies 
PLR,  Axy  =  1.1,  Y 

and 
=  3. 

Modal 

,5  .  . 

Loss 

Factors , 

>  •  •  •  •  o 

«  • 

• 

68 

19 

Modal  Frequencies 
PLR,  Axy  =  2.0,  Y 

and 
=  3. 

Modal 

,5  .  . 

Loss 

Factors , 

m  • 

ft 

69 

20 

Modal  Frequencies 
PLR,  Axy  =  4.0,  Y 

and 
=  3. 

Modal 

,5  .  . 

Loss 

•  •  1 

Factors, 

1  •  •  •  •  ^ 

9  9 

ft 

70 

21 

Modal  Frequencies 
PTR,  Axy  =  1.1,  Y 

and 
=  3, 

Modal 

.5  .  . 

Loss 

Factors , 

t  •  •  •  •  ^ 

o  o 

ft 

78 

22 

Modal  Frequencies 
PTR,  Axy  =  2.0,  Y 

and 
=  3, 

Modal 

.5  .  . 

Loss 

•  9  * 

Factors , 

ft  0 

ft 

79 

23 

Modal  Frequencies 
PTR,  Axy  =  4.0,  Y 

and 
=  3, 

Modal 

.5  .  . 

Loss 

Factors , 

ft  ft 

ft 

80 

24 

Modal  Frequencies 
PWU,  Axy  =  1.1,  Y 

and 
=  3, 

Modal 

.5  .  . 

Loss 

•  • 

Factors , 

•  ••••• 

.  e 

ft 

89 

25 

Modal  Frequencies 
PWU,  Axy  =  2.0,  Y 

and 
=  3, 

Modal 

.5  .  . 

Loss 

Factors , 

•  •  •  •  •  O 

•  • 

ft 

90 

26 

Modal  Frequencies 
PWU ,  Axy  =  4.0,  Y 

and 
=  3 

Modal 

.5  .  . 

Loss 

•  • 

Factors , 

»  •  •  •  •  o 

ft 

91 

27 

Modal  Frequencies 
PWR,  Axy  =  1.1,  Y 

and 
=  3 

Modal 

.5  «  . 

Loss 

•  « 

Factors , 

.  . 

»  ft 

99 

28 

Modal  Frequencies 
PWR ,  Axy  =  2.0,  Y 

and 
=  3 

Modal 

.5  .  . 

Loss 

Factors , 

•  •  •  •  •  9 

9  . 

9  ft 

100 

29 

Modal  Frequencies 
PWR,  Axy  =  4.0,  Y 

and 
=  3 

Modal 

.5  .  . 

Loss 

•  • 

Factors , 

•  ••••• 

,  O  , 

9  ft 

101 

30 

Results  from  Closed  Form  Solution  for 
Simply  Supported  Rectangular  Plate  .  .  .  . 

1  ft 

ft  ft 

121 

Aerospace  Structures 
Information  and  Analysis  Center 


1.0  INTRODUCTION 

The  use  of  viscoelastic  materials  for  vibration  control  has 
gained  wide  acceptance,  particularly  in  the  aerospace  industry. 
The  advantages  of  the  method  are  many  and  well  documented.  How¬ 
ever,  the  effectiveness  of  any  damping  treatment  is  critically 
dependent  on  its  geometry  and  on  the  properties  of  the  visco¬ 
elastic  material.  An  uninformed  choice  can  add  cost  and  weight 
but  fail  to  solve  the  vibration  problem.  Some  recent  advances  in 
analysis  techniques  utilize  finite  element  methods  to  provide  a 
more  reliable  and  systematic  approach  to  the  design  of  damping 
treatments  [1,2,3]* 

Plat  plate  sections  are  structural  elements  which  are  often 
good  candidates  for  integral  or  add-on  damping  treatments.  Such 
treatments  can  be  highly  effective  in  solving  problems  of 
fatigue,  acoustic  noise  radiation,  or  other  undesirable  effects 
of  resonant  vibration.  One  of  the  most  weight-efficient  forms  of 
damping  for  plate  applications  is  a  thin  layer  of  a  viscoelastic 
material  constrained  between  two  metal  face  sheets  to  form  a 
sandwich.  The  purpose  of  this  report  is  to  describe  several 
methods  for  the  analysis  and  design  of  sandwich  plate  structural 
elements. 

Three  methods  are  presented.  Each  represents  a  different 
trade-off  of  cost  of  use,  accuracy,  and  generality. 

The  first  and  most  general  technique  is  called  the  Modal 
Strain  Energy  (MSE)  method.  It  is  implemented  in  MSC/NASTRAN  and 
is  applicable  to  a  wide  range  of  structural  forms  in  addition  to 
plates.  Basically,  it  involves  modeling  the  viscoelastic  mater¬ 
ial  in  a  damped  structure  with  solid  elements  and  the  metallic 
material  with  solid,  plate,  shell,  or  other  elements  as  appropri¬ 
ate.  All  materials  are  treated  initially  as  being  purely  elastic 
(i.e.,  incapable  of  energy  dissipation).  Normal  mode  properties 
are  calculated  and  the  strain  energy  distribution  associated  with 
each  mode  shape  is  used  to  calculate  an  approximate  value  for  the 
modal  loss  factor.  The  theoretical  basis  of  the  MSE  method  is 
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described  later  in  this  report,  along  with  practical  considera¬ 
tions  for  modeling  of  sandwich  plates  using  MSC/NASTRAN.  The 
basic  assumptions  of  the  method  are  verified  by  comparisons  with 
an  exact  closed  form  solution  for  the  case  of  a  simply  supported, 
unriveted,  sandwich  rectangular  plate. 

The  second  method  is  essentially  a  set  of  design  charts  for 
sandwich  plates.  They  allow  a  designer  to  obtain  the  modal  fre¬ 
quencies  and  loss  factors  of  a  wide  variety  of  sandwich  plates 
with  only  simple  hand  calculations.  The  charts  were  compiled 
from  a  large  number  of  NASTRAN  runs  using  the  modal  strain  energy 
method.  They  are  plotted  in  dimensionless  form  for  generality 
and  give  modal  properties  for  the  first  four  modes  of  a  rectangu¬ 
lar  plate  with  various  boundary  conditions.  The  boundary  condi¬ 
tions  and  the  ranges  of  the  dimensionless  variables  are  chosen  to 
be  typical  of  the  situations  that  a  designer  might  commonly 
encounter  in  practice. 

The  third  method  is  a  simple,  inexpensive  technique  for 
designing  a  treatment  to  damp  the  higher  order  flexural  modes  of 
a  plate.  In  this  case  it  may  be  shown  that  the  boundary  condi¬ 
tions  are  relatively  unimportant.  A  closed  form  solution  is  used 
to  give  the  modal  loss  factor  as  a  function  of  natural  frequency. 
The  method  is  implemented  in  an  interactive  FORTRAN  program  and 
comparisons  are  made  with  NASTRAN  modal  strain  energy  results  to 
illustrate  the  effect  of  the  approximations. 
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2.0  MODAL  STRAIN  ENERGY  METHOD 

2.1  OVERVIEW 

The  discretized  equations  of  motion  for  a  damped  structure 
are  usually  written  in  the  form 

Mx  +  Cx  +  Kx  =  A(t)  ) 

where 

M,C,K  =  physical  coordinate  mass,  damping,  and  stiffness 
matrices  (all  real  and  constant) 

x,x,x  =  vectors  of  nodal  displacements,  velocities, 
and  accelerations 

Z  =  vector  of  applied  node  loads 

The  essence  of  the  modal  strain  energy  method  is  that  it  does  not 
attempt  to  find  the  damping  matrix  C.  This  would  he  impractical 
for  most  real  structures  and,  furthermore,  would  produce  a  system 
of  equations  which  would  he  very  costly  to  solve.  Rather,  in  the 
MSE  method,  one  assumes  that  the  damped  structure  can  he  repre¬ 
sented  in  terms  of  the  real  normal  modes  of  the  associated 
undamped  system  if  appropriate  damping  terms  are  inserted  into 
the  uncoupled  modal  equations  of  motion.  That  is: 

“j.  +  =  ^j.(t)  (2) 

r  =  1,2,5  (5) 

~  r 


where 

Oj.  =  r'th  modal  coordinate 

Uj.  =  natural  radian  frequency  of  the  r'th  mode 

^(r)  _  pf-th  mode  shape  vector  of  the  associated  undamped 
system 

n(r)  =  loss  factor  of  the  r'th  mode 
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Equations  (2)  and  (3)  imply  that  the  physical  coordinate  damping 
matrix  C  of  Eq.  (1)  need  not  be  explicitly  calculated  but  that  it 
can  be  diagonalized,  at  least  approximately,  by  the  same  real 

modal  matrix  that  diagonalizes  K  and  M. 

The  modal  loss  factors  are  calculated  by  using  the  undamped 

mode  shapes  and  the  material  loss  factor  for  each  material.  For 
structures  containing  a  viscoelastic  material,  the  material  loss 
factor  of  the  metal  is  very  small  compared  to  that  of  the  visco¬ 
elastic.  In  this  situation  the  modal  loss  factor  is  found  from 


where 


material  loss  factor  of  viscoelastic  evaluated 
at  the  r’th  calculated  resonant  frequency 


fraction  of  elastic  strain  energy  attributable 
to  the  viscoelastic  when  the  structure  deforms 
in  the  r'th  mode  shape 


A  derivation  of  Eq.  (4)  is  given  in  Section  2.2  of  this 
report.  It  is  shown  that  modal  loss  factors  obtained  from 
Eq.  (4)  will  approximate  those  obtained  from  the  complex 
stiffness  eigenvalues  of  the  complementary  solution  of  Eq.  (l). 
However,  the  modal  strain  energy  approach  has  the  advantage  of 

much  lower  cost. 

Calculation  of  the  modal  energy  distributions  fits  quite 
naturally  within  finite  element  methods  and  is  a  standard  option 
in  some  commercial  codes.  Further  advantages  of  the  method  are 
that  only  undamped  normal  modes  need  be  calculated  and  that  the 
energy  distributions  obtained  are  of  direct  use  to  the  designer 
in  deciding  where  to  locate  damping  material.  The  disadvantage 
is  that  some  approximation  is  required  to  accommodate  the 
frequency-dependent  properties  commonly  found  in  viscoelastic 

materials . 
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2 . 2  THEORY 

An  approximate  expression  is  derived  below  for  the  modal 
loss  factor  obtained  from  an  eigenvalue  analysis  of  a  structure 
with  complex  stiffness. 

One  form  of  the  discretized  (i.e.,  finite  element)  version 
of  a  partial  differential  equation  for  free  vibration  of  a 
structure  is: 

Mx+Kx=0  (5) 


where  the  stiffness  matrix  K  is  constant  but  complex  if  the 
structure  contains  a  viscoelastic  material.  Equation  (5)  is 
converted  to  an  eigenvalue  problem  by  assuming  a  solution  of  the 
form 


X  = 


ip? 

e 


(6) 


where  p*  and  are  the  r'th  complex  eigenvalue  and  eigen¬ 

vector.  That  is. 


=  4^)  +  i4-)  (7 

P?  = 

where  and  p^  are  real.  The  term  is  the 

loss  factor  for  the  r'th  mode.  The  eigenvalue  problem  is  then, 
from  Eqs.  (5)  and  (6): 


K  r  =  P*  ^ 


(9) 


and  p*  would  be  real  and 
■^r 


Now  if  K  were  purely  real,  ^ 
related  by  the  usual  Rayleigh's  quotient  formula: 


rrT 


(10) 
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where  the  *  superscript  is  dropped  to  denote  a  real  quantity. 

If  K  is  perturbed  by  6K,  where  6K  is  complex,  Pj.^  will  likewise 

^  2 
acquire  an  imaginary  part  which  may  be  written  as  inp  after 

Eq.  (8).  Then,  if  the  perturbed  stiffness  matrix  is  written  as 


K  =  K„  +  i  K, 


(11) 


the  following  is  obtained  from  Eqs.  (8),  (10),  and  (11),  after 
dropping  the  mode  index  r 


p^d+in) 


K  $>* 

- m S -  +  i  - m - 

M  M 


(12) 


An  approximate  value  for  n  can  be  calculated  by  approximating  the 
complex  eigenvector  by  the  real  vector  calculated  from 
purely  elastic  analysis,  i.e.,  by  suppressing  the  imaginary  part 
of  K  .  The  approach  is  essentially  an  extension  of  Rayleigh’s 
principle  into  the  complex  domain.  Making  this  approximation  in 
Eq.  (12)  and  equating  real  and  imaginary  parts  gives 

p2  ,  (13) 

M  $ 


P 


“T? 

$  M  $ 

rsj 


(14) 


If  the  matrix  K  is  obtained  by  finite  element  analysis,  it 
may  be  divided  into  two  additive  terms.  The  first,  called  K  ,  is 
obtained  from  contributions  of  the  purely  elastic  elements  (the 
metallic  portion  of  the  structure).  The  second,  called  K^,  is 
obtained  from  the  solid  elements  (used  to  model  the  viscoelastic 
material).  Both  terms  are  matrices  of  the  same  order  as  K, 


K 


K 


K 


-V 


(15) 
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K  will  "be  completely  real.  will  be  complex  but,  for  the 

present  case  where  only  a  single  viscoelastic  material  is 
involved,  its  imaginary  and  real  parts  will  have  the  ratio 
where  is  the  material  loss  factor  of  the  core.  Then, 


=  SvR  +  1  Svi 


=  EvR 


(16) 

(17) 


By  previous  assumption,  only  contributes  to  Kj  so 
=  K 


-vl 


(18) 


When  a  purely  elastic  normal  modes  analysis  is  performed,  the 
strain  energy  associated  with  a  given  mode  shape  is 


V  =  1  J 


(19) 


The  portion  of  this  energy  which  is  attributable  to  strain  in  the 
core  is 


Eliminating  p^  between  Eq.  (12)  and  (13)  gives 


(20) 


n  =  n 


Si  i. 
''  i’’  Kr  t 


(21) 


Combining  Eq.  (17)  through  (21)  and  reinstating  the  mode  index 
superscript  gives  the  final  result  for  modal  loss  factor  in  terms 
of  elastic  energies 


(r)  n  V 

\  TTFT 


(22) 
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This  derivation  is  intended  to  motivate  and  clarify  the 
comparison  of  results  from  complex  eigenvalue  analysis  and  modal 
strain  energy  analysis.  It  should  he  noted,  however,  that  the 
problem  statement  itself  is  not  entirely  realistic.  It  is  well 
known  that  complex  stiffnesses  that  do  not  vary  with  frequency 
lead  to  system  responses  that  are  noncausal  (response  anticipates 
input)  and  hence  are  not  physically  realizable.  Nonetheless,  the 
comparison  is  believed  to  be  useful  in  that  the  symptoms  of 
noncausality  are  quite  weak  [4]  and  the  identical  assumptions  are 
applied  in  both  methods. 


2.5  FINITE  ELEMENT  MODELING  OP  THREE-LAYER  PLATES 
2.3.1  Choice  of  Elements 

Modeling  of  sandwich  structures  requires  that  the  strain 
energy  due  to  shearing  of  the  core  be  accurately  represented. 
Practical  considerations  dictate  that  this  be  done  with  minimum 
increase  in  computation  cost  relative  to  a  uniform,  single-layer 
model.  In  this  section,  a  modeling  method  is  described  which  is 
reasonably  efficient  and  has  the  important  advantage  of  being 
readily  implemented  within  MSC/NASTRAN,  a  widely  available  code. 

Figure  1  shows  the  arrangement  for  modeling  of  a  three-layer 
sandwich.  The  face  sheets  are  modeled  with  quadrilateral  or 
triangular  plate  elements  producing  stiffness  at  two  rotational 
and  three  translational  degrees  of  freedom  per  node.  The  visco¬ 
elastic  core  is  modeled  with  solid  elements  producing  stiffness 
at  three  translational  degrees  of  freedom  per  node.  All  nodes 
are  at  element  corners.  In  MSC/NASTRAN,  the  plate  elements  are 
called  TRIA3's,  QUAD4*s,  TRIA6's,  and  QUADS 's,  and  the  solid 
elements  are  called  PENTA’s  and  HEXA's.  A  key  feature  of  these 
plate  elements  in  the  present  application  is  their  ability  to 
account  for  coupling  between  stretching  and  bending  deformations 
[5].  This  allows  the  plate  nodes  to  be  offset  to  one  surface  of 
the  plate,  coincident  with  the  corner  nodes  of  the  adjoining 
solid  elements.  In  this  way  a  three-layer  plate  can  be  modeled 
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Figure  1 


{  UPPER  FACE  j 
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3 
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LOWER  FACE  \ 

I'l 

SHEET  1  , 

PANEL  SIDE  VIEW 

FINITE  ELEMENT 
REPRESENTATION 


Finite  element  modeling  of  a  sandwich 
panel  with  viscoelastic  core 
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with  only  two  layers  of  nodes.  Earlier  methods  implemented 
within  NASTRAN  were  restricted  to  beams  [6]  or  required  four 
layers  of  nodes  and  extensive  constraint  equations  to  achieve  the 
proper  bending-shearing  behavior  of  the  sandwich  [7].  Aspect 
ratios  of  the  solid  elements  (in-plane  dimension/thickness 
dimension)  as  high  as  5000  have  been  used  successfully  to  model 
the  thin  viscoelastic  core  layers.  In  all  analyses  reported 
here,  Poisson’s  ratio  of  the  core  elements  is  taken  to  be  0.49* 

2.5.2  Reduction  of  Equations  of  Motion 

In  all  but  the  smallest  problems,  the  mass  and  stiffness 
matrices  are  condensed  by  partitioning  and  performing  a  Guyan 
reduction  prior  to  calculation  of  eigenvalues.  As  usual  in 
vibration  analysis,  some  care  is  warranted  in  the  selection  of 
the  degrees  of  freedom  to  be  retained  during  this  reduction.  Eor 
the  sandwich  structures  analyzed  in  this  report,  only  out-of- 
plane  displacements  need  be  retained.  Some  displacements  should 
be  retained  for  both  face  sheets,  although  it  is  not  necessary  to 
jfe0p  "both  upper  and  lower  face  displacements  at  any  single  loca¬ 
tion  on  the  model.  If  out-of-plane  displacements  of  only  one 
face  sheet  are  kept,  the  results  for  natural  frequencies  as  well 
as  core-to-total  energy  ratios  can  show  a  pronounced  dependence 
on  the  Poisson's  ratio  of  the  core.  Although  such  a  dependence 
is  probably  real  for  some  cases,  such  as  doubly  curved  shells,  it 
should  not  occur  for  simpler  cases  such  as  straight  sandwich 
■beams — and  in  fact  does  not  occur  if  the  rule  given  above  is 
observed  in  reducing  the  discretized  equations  of  motion.  Exist¬ 
ing  data  on  Poisson's  ratio  of  most  viscoelastic  materials  are 
probably  not  adequate  for  accurate  modeling  of  doubly  curved 
sandwich  shells  in  the  important  transition  region  of  the 
material,  and  certainly  not  in  the  glassy  region. 
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2.3*3  Solution  Method 

In  the  modal  strain  energy  method,  a  standard  normal  mode 
extraction  run  is  made  with  all  material  constants  treated  as 
real  and  constant.  The  elastic  strain  energy  in  each  element  for 
each  mode  is  calculated,  as  well  as  the  energy  fraction  in  the 
viscoelastic  core  for  each  mode.  These  fractions  multiplied  by 
the  core  material  loss  factor  give  the  modal  loss  factors  which 
are  input  via  a  damping  vs.  frequency  table  for  use  in  subsequent 
forced  response  calculations. 

A  basic  difficulty  with  the  modal  strain  energy  method  (or 
any  normal  mode  method)  is  that  the  modal  properties  are  obtained 
from  system  matrices  that  are  assumed  to  be  constant.  Visco¬ 
elastic  materials,  however,  have  storage  moduli  which  vary  signi¬ 
ficantly  with  frequency.  There  is  no  theoretically  correct  way 
to  resolve  this  contradiction.  However,  there  are  great  practi¬ 
cal  advantages  to  making  response  predictions  in  terms  of  a 
normal  mode  set  obtained  from  constant  material  properties.  This 
can  be  done  with  reasonable  accuracy  if  a  simple  correction  is 
made  to  the  modal  loss  factors  obtained  by  Equation  (4).  This 
correction  is  only  to  the  modal  damping  ratios  because  these  are 
the  only  modal  parameters  that  can  be  readily  adjusted  by  the 
finite  element  analyst.  It  is  explained  here  for  completeness 
even  though  no  forced  response  calculations  were  performed  for 
this  report.  The  correction  is  obtained  as  follows. 

For  broadband  excitation,  most  of  the  response  of  a  given 
mode  occurs  within  a  narrow  band  around  the  mode's  natural  fre¬ 
quency.  It  is  natural  then  to  require  that  the  energy  distribu¬ 
tion  used  to  compute  the  loss  factor  for  a  given  mode  be  obtained 
using  a  stiffness  matrix  evaluated  for  material  properties  taken 
at  that  mode's  frequency.  Because  the  natural  frequencies  them¬ 
selves  depend  on  material  properties,  an  iterative  solution  of 
two  simultaneous  relations  (the  eigenvalue  problem  for  each  mode 
number  and  the  material  property  vs.  frequency  relation)  is  re¬ 
quired.  This  is  readily  done  [1],  but  a  further  problem  remains. 
The  final  modal  coordinate  representation  of  the  structure  must 
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come  from  a  single  stiffness  matrix  evaliiated  using  a  single 
value  of  storage  modulus  for  the  core  material.  Natural  fre¬ 
quencies,  mode  shapes,  and  modal  masses  will  he  correct  for,  at 
most,  one  mode.  A  further  correction  of  the  modal  loss  factor 
has  heen  found  to  give  some  improvement. 

Each  modal  equation  of  motion  has  the  form  given  in  Equation 
(2).  At  resonance  the  first  and  last  terms  on  the  left  cancel 
each  other.  The  response  magnitude  is  inversely  proportional  to 
the  product  which  is  the  coefficient  of  a^,  the  modal 

velocity.  If  is  altered  to  correct  for  the  error  in  an 

improvement  in  peak  response  can  he  expected,  although  resonance 
will  still  occur  at  a  slightly  shifted  frequency  and  some  error 
will  remain  due  to  which  depends  on  modal  mass.  In  test  cases 
run  for  sandwich  beams  [1],  it  was  found  that  taking  to  he 
proportional  to  \/(t2  (Gr2  =  core  shear  modulus)  would  improve  the 
agreement  between  the  MSE  method  and  the  direct  frequency 
response  method.  This  is  of  course  an  approximation  since 
depends  on  properties  of  the  face  sheets  as  well  as  the  core. 

The  modal  damping  ratios  are  adjusted  according  to 


n 


(r)' 


/ 


M£rl 

^2, ref 


(25) 


where 

=  adjusted  modal  damping  ratio  for  the  r^th  mode 

=  modal  damping  ratio  for  the  r'th  mode  obtained 
by  iteration 


®2,ref 


core  shear  modulus  used  in  final  normal  modes 
calculation  to  obtain  modal  frequencies,  shapes, 
and  masses 

core  shear  modulus  at  f  =  f^  where  f„  is  r'th 

r  r 

mode  frequency  calculated  with  =  G2 
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2.4  EXAMPLE 

A  closed  form  solution  exists  for  the  complex  eigenvalues 
(i.e.,  natural  frequencies  and  modal  loss  factors)  of  a  simply 
supported  sandwich  plate  [8].  The  solution  is  described  in 
Section  4.0  of  this  report  in  connection  with  the  design  of  con¬ 
strained-layer  damping  treatments  for  high-order  local  modes  of 
plate  sections.  In  this  section  it  is  used  to  verify  one  of  the 
most  important  implications  of  the  MSE  method;  namely,  that  the 
modal  loss  factors  for  all  modes  of  a  sandwich  plate  are  directly 
proportional  to  the  material  loss  factor  of  the  viscoelastic 
core.  The  sample  problem  is  also  used  to  illustrate  the  input 
data  to  NASTRAN  for  MSE  analysis  of  a  sandwich  plate. 

Figures  2  through  5  give  a  comparison  of  modal  loss  factors 
obtained  by  using  MSC/NASTRAN  and  MSE  (MSC/NASTRAN-MSE)  and  by 
the  closed  form  solution  of  Ref.  [8].  The  format  of  these  plots 
is  used  throughout  this  report  and  is  explained  in  detail  in 
Section  3.1.  In  brief,  the  ordinate  is  a  dimensionless  quantity 
proportional  to  modal  loss  factor  and  the  abscissa  is  a  dimen¬ 
sionless  quantity  proportional  to  the  shear  modulus  of  the  visco¬ 
elastic  core.  The  curves  marked  with  specific  values  of  are 
obtained  from  the  closed  form  solution.  The  remaining  curve, 
obtained  by  the  MSC/NASTRAN-MSE  method,  gives  results  that  are 

inherently  independent  of  n^. 

The  test  case  is  a  simply  supported  rectangular  sandwich 

plate  of  the  following  dimensions: 


in-plane  dimensions 
upper  face  sheet  thickness 
lower  face  sheet  thickness 
core  thickness 
face  sheet  material 


=  10”  X  11" 

=  0.055" 

=  0.055" 

=  0.0045" 

=  aluminum 
E  =  1  o"^  psi 
P  =  0.1  Ib/in^ 
V  =  0.3 


shear  modulus  of  core  material  =  variable 
loss  factor  of  the  core  material  =  variable 
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Figures  2  through  3  show  that  the  closed  form  and 
MSC/NASTRAN-MSE  results  agree  very  closely  for  small  values  of 
the  material  loss  factor.  Some  divergence  is  seen  for  larger 
values  on  the  order  of  unity  or  greater.  The  agreement  also 
depends  on  the  value  of  the  shear  parameter  g.  It  is  best  for  g 
equal  to  or  less  than  the  value  giving  highest  damping. 
Fortunately,  most  practical  constrained  layer  treatments  tend  to 
fall  in  this  range. 

A  tabular  representation  of  the  closed  form  results  used  to 
prepare  Figures  2  through  3  is  given  in  Tables  1  through  4* 
Results  for  higher  modes  are  also  given  in  the  tables. 

Sample  NASTRAN  input  and  output  are  given  in  Appendix  A  for 
a  plate  with  the  properties  listed  above  and  a  core  shear  modulus 
of  450  psi.  This  sample  case  corresponds  to  a  dimensionless 
shear  parameter  of  g  =  40. 
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Figure  2  Damping  of  the  first  mode  of  simply  supported  rectangular  sandwich 
plate  obtained  by  NASTRAN/Modal  Strain  Energy  method  and  by 
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Figure  3  Damping  of  the  second  mode  of  simply  sunported  rectangular  sandwich 
plate  obtained  by  NASTRAN/Modal  Strain  Energy  method  and  by 
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Figure  4  Damping  of  the  third  mode  of  simply  supported  rectangular  sandwich 
plate  obtained  by  NASTRAN/Modal  Strain  Energy  method  and  by 
exact  complex  eigenvalue  solution  [8] 
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Figure  5  Damping  of  "the  fourth  mode  of  simply  supported  rectangular  sandwich 
plate  obtained  by  NASTRAN/Modal  Strain  Energy  method  and  by 
exact  complex  eigenvalue  solution  [8] 
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3.0  DESIGN  CHARTS  FOR  SANDWICH  PLATES 

In  this  section,  sets  of  design  charts  are  given  which  allow 
a  user  to  rapidly  estimate  the  loss  factors  and  natural 
frequencies  for  a  rectangular  sandwich  plate  of  various  boundary 
conditions.  The  charts  were  compiled  from  a  large  number  of 
NASTRAN  analyses  using  the  modal  strain  energy  method.  They 
allow  fairly  accurate  predictions  of  damping  which  take  boundary 
conditions  into  account  and  yet  do  not  require  the  user  to  actu¬ 
ally  prepare  or  run  any  finite  element  models.  The  usefulness  of 
these  charts  derives  from  the  fact  that,  to  the  authors’  know¬ 
ledge,  no  exact  solutions  exist  for  other  than  simply  supported 
boundary  conditions. 


3.1  DATA  FORMAT 

The  charts  are  in  terms  of  dimensionless  variables  in  order 
to  convey  the  maximum  amount  of  information.  It  may  be  shown  [9] 
that  a  rectangular  sandwich  plate  can  be  completely  described  by 
four  dimensionless  parameters: 


=  core  material  loss  factor 


g  =  shear  parameter 

Tt  /  1  _  1  ^ 

""  T^  ^  Vi  ^  ) 


Y  =  geometry  parameter 

(T,+Tj+2T2)^  ^ 

"  4D(1-v") 

Axy  =  in-plane  aspect  ratio 
=  b/a 


(24) 


(25) 

(26) 
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where 


T,,!, 


E,  ,E 

a,b 


3 


D 

V 


thicknesses  of  the  face  sheets 

thickness  of  the  core  layer 

real  part  of  the  complex  shear  modulus 
[Tj(1+in^)]  of  the  viscoelastic  material 

Young's  moduli  of  the  face  sheets 

in-plane  dimensions  of  the  plate 

sum  of  the  flexural  stiffnesses  of  the  upper  and 
lower  face  sheets,  each  about  its  own  center  plane 

Poisson’s  ratio  of  the  face  sheets 


(  I*  )  t 

The  charts  give  the  dimensionless  loss  parameter,  n  /\» 
as  a  function  of  the  shear  parameter  for  the  first  four  bending 
modes  of  a  plate  for  various  values  of  the  aspect  ratio  and  geo¬ 
metry  parameter.  A  value  of  approximately  Y  =  3*5  characterizes 
a  sandwich  plate  with  equal  thickness  face  sheets.  The  situation 
of  equal  face  sheets  is  fairly  common  in  practice  and  therefore 
is  included  for  all  boundary  conditions.  Additional  values  of  Y 
are  included  for  some  boundary  conditions  to  cover  the  case  of 
unequal  face  sheets  which  often  occurs  with  add-on  constrained 
layer  damping  treatments. 

Natural  frequencies,  in  a  normalized  form,  are  also  given  as 
a  function  of  shear  parameter  for  the  first  four  modes.  The  form 
is  f_/fQ.j  where  f^.  is  the  natural  frequency  of  the  r'th  mode. 

The  reference  frequency  fQ^  is  defined  as  the  first  natural  fre¬ 
quency  of  a  simply  supported  plate  of  the  same  in-plane  dimen¬ 
sions  as  the  actual  plate  but  with  flexural  stiffness  equal  to 
the  sum  of  the  stiffnesses  of  the  upper  and  lower  face  sheets. 
They  are  calculated  using  the  formula  [9]: 

where 

p  =  mass  density  per  unit  area  of  the  plate 
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The  reference  frequencies  calculated  by  Eq.  (27)  are  given 
in  Table  5.  It  may  be  noted  that  they  are  not  strictly  functions 
of  only  the  shear  and  geometry  parameters  but  also  depend  on  the 
material  properties  of  the  face  sheets.  Since  two  different  sets 
of  properties  (corresponding  to  steel  and  aluminum)  were  used  to 
span  the  desired  range  of  the  shear  parameter,  the  reference 
frequencies  corresponding  to  each  are  given. 

Summary  tables  of  the  NASTRAN  results  used  to  prepare  the 
frequency  and  damping  plots  are  given  for  each  set  of  boundary 
conditions.  The  tables  give  results  for  the  fifth  and  higher 
modes,  in  addition  to  the  first  four  modes  for  which  results  are 
plotted. 

Each  graph  and  table  is  marked  with  a  three  letter  abbrevia¬ 
tion  denoting  the  boundary  conditions.  The  first  character  in 
every  case  is  a  P  (pinned)  and  indicates  that  out— of— plane  dis¬ 
placements  of  both  face  sheets  were  constrained.  The  second 
character  is  either  a  T  (tilting),  L  (level),  or  V  (wind-up). 
Respectively,  they  designate  an  unconstrained,  perfectly  con¬ 
strained,  or  elastically  constrained  condition  on  rotation  of  the 
face  sheets  about  an  axis  parallel  to  the  plate  edge.  The  third 
character  is  either  a  U  (unriveted)  or  an  R  (riveted).  Riveted 
implies  that  shearing  deformation  of  the  core  has  been 
constrained  along  the  plate  edge. 
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for  shear  parameter  =  1000. 
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3.2  DESIGN  CHARTS 

3.2.1  PTU  Boundary  Conditions 

The  results  for  PTU  (simply-supported,  unriveted)  boundary 
conditions  are  given  in  Figures  6  through  14  and  Tables  6  through 
17.  These  boundary  conditions  are  likely  to  be  appropriate  for 
plate  sections  in  lightweight,  built-up  structures  using  add-on 
constrained  layer  damping.  Several  values  of  the  geometry 
parameter  are  used  since  it  is  likely  to  be  a  design  variable  in 
these*  situations. 

Figure  6  shows  damping  as  a  function  of  shear  parameter  for 
the  first  four  modes,  with  a  geometry  parameter  of  Y  =  3*5  and  an 
aspect  ratio  Axy  =  1.1.  Figures  7  and  8  show  similar  information 
for  aspect  ratios  of  2.0  and  4.0.  The  next  six  figures,  9 
through  14,  give  similar  data  but  organized  to  show  how  damping 
varies  with  the  geometric  parameter  as  well  as  the  shear  para¬ 
meter.  For  clarity,  the  data  covering  four  modes  (for  each  value 
of  aspect  ratio)  is  split  into  two  plots,  with  each  plot  covering 
two  alternate  modes. 

Natural  frequencies  of  sandwich  plates  with  PTU  boundary 
conditions  can  be  obtained  from  Figures  15  through  26.  Each  plot 
gives  results  for  one  of  three  aspect  ratios  (1.1,  2.0,  or  4.0) 
and  one  of  four  geometry  parameters  (0.5,  1.5,  5.5,  or  4.5). 
Reference  frequencies  used  in  normalizing  the  data  of  these 
figures  are  given  in  Table  5. 

A  tabular  representation  of  the  data  in  Figures  6  through  26 
as  well  as  data  for  higher  modes  is  given  in  Tables  6  through  17. 
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Figure  6  Damping  of  a  sandwich  rectangular  plate,  PTU  boundary  conditions, 
Axy  =  1.1,  Y  =  3.5 
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Figure  8  Damping  of  a  sandwich  rectangular  plate,  PTU  boundary  conditions, 
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Figure  9  Damping  of  a  sandwich  rectangular  plate,  PTU  boundary  conditions, 
Axy  =  1.1,  modes  1  and  3,  variable  Y  and  g 
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Figure  10  Damping  of  a  sandwich  rectangular  plate,  PTU  boundary  conditions 
Axy  =  1.1,  modes  2  and  4,  variable  Y  and  g 
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Figure  11  Damping  of  a  sandwich  rectangular  plate,  PTU  boundary  conditions, 
Axy  =  2.0,  modes  1  and  3,  variable  Y  and  g 
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Figure  12  Damping  of  a  sandwich  rectangular  plate,  PTU  boundary  conditions 
Axy  =  2.0,  modes  2  and  4,  variable  Y  and  g 
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Figure  13  Damping  of  a  sandwich  rectangular  plate,  PTU  boundary  conditions 
Axy  =  4.0,  modes  1  and  3,  variable  Y  and  g 
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Figure  14  Damping  of  a  sandwich  ■rectangular  plafe,  PTU  boundary  conditions, 
Axy  =  4.0,  modes  2  and  4,  variable  Y  and  g 


Aerospace  Structures 


Information  and  Analysis  Center 


CO 

o 


CN 

o 


tH 

O 


O 

O 

tH 


bO 


U 

CD 

-P 

CD 

B 

p 

fd 

Pu 

u 

fd 

CD 

rC 

CO 


a9q.auipapd  Aouanbaaj 


37 


Figure  15  Natural  frequencies  of  sandwich  rectangular  plate 
PTU  boundary  conditions,  Axy  =  1.1,  Y  =0.5 
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Figure  16  Natural  frequencies  of  sandwich  rectangular  plate, 
PTU  boundary  condition,  Axy  =  2.0,  Y  =  0.5 
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Figure  17  Natural  frequencies  of  sandwich  rectangular  plate 
PTU  boundary  conditions ,  Axy  =  4.0,  Y  =  0.5 
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Figure  18  Natural  frequencies  of  sandwich  rectangular  plate 
PTU  boundary  conditions,  Axy  =  1.1,  Y  =  1.5 
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Figure  19  Natural  frequencies  of  sandwich  rectangular  plate 
PTU  boundary  conditions ,  Axy  =  2 . 0 ,  Y  =  1 . 5 
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Figure  20  Natural  frequencies  of  sandwich  rectangular  plate, 
PTU  boundary  conditions,  Axy  =  4.0,  Y  =1.5 
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Figure  21  Natural  frequencies  of  sandwich  plate,  PTU 
boundary  conditions-,  Axy  =  1.1,  Y  =  3.5 
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Figure  22  Natural  frequencies  of  sandwich  rectangular  plate 
PTU  boundary  conditions ,  Axy  =  2.0,  Y  =  3.5 
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Figure  23  Natural  frequencies  of  sandwich  rectangular  plate, 
PTU  boundary  conditions,  Axy  =  4.0,  Y  =  3.5 


Aerospace  Structures 
Information  and  Analysis  Center 


o  o  o  o 

•  •  •  • 

00  VD  ^  CVJ 


CO 

o 


bO 


o 


U 

-M 

0) 

B 

fd 


(d 

0) 

rC 

CO 


tH 

o 

tH 


aaq.auieaPd  Aouanbaaj 
46 


0) 

-H 

fd 

iH 

Ph 

P  LO 
fd  .• 

rH 

P 

bO  II 

C 

fd 

•H 

O 

(U  vH 

P  • 

vH 

O  II 
•H 

15 
T3 
P 
fd 
CO 

4^ 

O 

CO 

<u 

•H 

a 
P 
(U 
P 
cr* 

<U 

4h 

I— I 
fd 
P 
P 
-p 
fd 


ZJ- 

CM 

cu 

p 

p 

bO 

•H 

Pm 


PTU  boundary  conditions,  Axy 


Aerospace  Structures 


asq-auipapj  Aouanbaaj 
47 


Figure  25  Natural  frequencies  of  sandwich  rectangular  plate, 
PTU  boundary  conditions,  Axy  =  2.0,  Y  =4.5 
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Figure  26  Natural  frequencies  of  sandwich  rectangular  plate 
PTU  boundary  conditions ,  Axy  =  4.0,  Y  =  4.5 


TABLE  6 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PTU  (zero  translation,  unrestrained  rotation,  unrestrained  shear) 
Aspect  Ratio  (Axy)  =  !•! 

Geometric  Parameter  (Y)  =  0.5 
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TABLE  7 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PTU  (zero  translation,  unrestrained  rotation,  unrestrained  shear) 
Aspect  Ratio  (Axy)  =  1,1 
Geometric  Parameter  (Y)  =  1.5 
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7167.  0.009  4258.  0.0455  748.  0.120  884.  0.191  993.  0.090 


TABLE  10 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PTU  (zero  translation,  unrestrained  rotation,  unrestrained  shear) 
Aspect  Ratio  (Axy)  =  2.0 
Geometric  Parameter  (Y)  =0.5 
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TABLE  13 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PTU  (zero  translation,  unrestrained  rotation,  unrestrained  shear) 
Aspect  Ratio  (Axy)  =2.0 
Geometric  Parameter  (Y)  =4.5 


Aerospace  Structures 
Information  and  Analysis  Center 


0.212 


TABLE  14 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PTU  (zero  translation,  unrestrained  rotation,  unrestrained  shear) 
Aspect  Ratio  (Axy)  =  4.0 
Geometric  Parameter  (Y)  =  0.5 
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MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PTU  (zero  translation,  unrestrained  rotation,  unrestrained  shear) 
Aspect  Ratio  (Axy)  =  4.0 
Geometric  Parameter  (Y)  =3.5 
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3.2.2  PLR  Boundary  Conditions 

Damping  properties  of  sandwich  plates  with  PLR  (fixed) 
boundary  conditions  are  given  in  Figures  27  through  29*  Only  one 
value  of  the  geometry  parameter,  Y  =  3»5»  is  considered.  This 
corresponds  to  the  situation  of  equal  face  sheet  thicknesses. 

Natural  frequencies  of  sandwich  plates  with  PLR  boundary 
conditions  are  given  in  Figures  30  through  32.  Reference 
frequencies  are  given  in  Table  3* 

A  tabular  presentation  of  the  data  for  damping  and  natural 
frequencies  of  plates  with  PLR  boundary  conditions  is  given  in 
Tables  18  through  20. 
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Figure  27  Damping  of  a  sandwich  rectangular  plate 
PLR  boundary  conditions,  Axy  =  1.1,  Y  = 
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Figure  28  Damping  of  a  sandwich  rectangular  plate, 
PLR  boundary  conditions,  Axy  =  2.0,  Y  = 
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Figure  29  Damping  of  a  sandwich  rectangular  plate, 
PLR  boundary  conditions ,  Axy  =  4.0,  Y  = 
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Figure  30  Natural  frequencies  of  a  rectangular  sandwich 
plate,  PLR  boundary  conditions,  Axy  =  1.1,  Y 
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Figure  31  Natural  frequencies  of  a  rectangular  sandwich 
plate,  PLR  boundary  conditions,  Axy  =  2.0,  Y 
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Figure  32  Natural  frequencies  of  a  rectangular  sandwich 
plate,  PLR  boundary  conditions,  Axy  =  4.0,  Y 


TABLE  18 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PLR  (zero  translation,  zero  rotation,  zero  shear) 
Aspect  Ratio  (Axy)  =1.1 
Geometric  Parameter  (Y)  =  3.5 
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0.092  1269.  0.205  1451.  0.271  1908.  0.152 


TABLE  19 

MODAL  FREQUENCIES  AND  MODAL  LOSS. FACTORS 

Boundary  Condition  =  PLR  (zero  translation,  zero  rotation,  zero  shear) 
Aspect  Ratio  (Axy)  =  2.0 
Geometric  Parameter  (Y)  =  3.5 
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TABLE  20 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PLR  (zero  translation,  zero  rotation,  zero  shear) 
Aspect  Ratio  (Axy)  =  4.0 
Geometric  Parameter  (Y)  =  3.5 
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5663.  0.174  6999.  0.259 


Aerospace  Structures 
Information  and  Analysis  Center 


3.2.3  PTR  Boundary  Conditions 

Damping  properties  of  sandwich  plates  with  PTR  (simply- 
supported,  riveted)  boundary  conditions  are  given  in  Figures  33 
through  35.  Again,  only  a  single  value  of  the  geometry 
parameter,  Y  =  3*5>  is  considered. 

Natural  frequencies  of  sandwich  plates  with  PTR  boundary 
conditions  are  given  in  Figures  36  through  Reference 

frequencies  are  given  in  Table  5* 

A  tabular  presentation  of  the  data  for  damping  and  natural 
frequencies  of  plates  with  PTR  boundary  conditions  is  given  in 
Tables  21  through  23* 
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Figure  33  Damping  of  a  sandwich  rectangular  plate,  PTR 
boundary  conditions,  Axy  =  1.1,  Y  =  3.5 
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Figure  34  Damping  of  a  sandwich  rectangular  plate,  PTR 
boundary  conditions,  Axy  =2.0,  Y  =  3. '5 
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Figure  35  Damping  of  a  sandwich  rectangular  plate,  PTR 
boundary  conditions ,  Axy  =  4.0,  Y  =  3.5 
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Figure  36  Natural  frequencies  of  a  rectangular  sandwich 
plate ,  PTR  boundary  conditions ,  Axy  =  1.1,  Y 
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Figure  37  Natural  frequencies  of  a  rectangular  sandwich 
plate,  PTR  boundary  conditions ,  Axy  =  2.0,  Y 
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Figure  38  Natural  frequencies  of  a  rectangular  sandwich 
plate,  PTR  boundary  conditions,  Axy  =  4.0,  Y 


TABLE  21 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PTR  (zero  translation,  unrestrained  rotation,  zero  shear) 
Aspect  Ratio  (Axy)  =  1.1 
Geometric  Parameter  (Y)  =  3.5 
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TABLE  22 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PTR  (zero  translation,  unrestrained  rotation,  zero  shear) 
Aspect  Ratio  (Axy)  =  2.0 
Geometric  Parameter  (Y)  =  3.5 
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TABLE  23 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PTR  (zero  translation,  unrestrained  rotation,  zero  shear) 
Aspect  Ratio  (Axy)  =  4.0 
Geometric  Parameter  (Y)  =  3.5 
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5.2.4  PWU  Boundary  Conditions 

A  situation  that  might  lead  to  PV/U  or  PWR  boundary  condi¬ 
tions  is  shown  in  Figure  39*  A  structure  is  fabricated  by  butt 
welding  of  plate  sections  that  contain  an  integral  damping  treat¬ 
ment.  One  leg  of  the  weldment  sees  a  constraint  on  shearing  of 
the  sandwich  core  (PWR),  while  the  other  [PWU]  does  not.  Both 
see  some  restraint  on  bending  rotation  at  the  welded  boundary  but 
it  is  not  held  exactly  to  zero.  The  restraint  is  approximated  as 
a  pure  stiffness  and  is  evaluated,  somewhat  arbitrarily,  as 
follows.  The  degree  of  elastic  restraint  is  taken  to  be  equal  to 
the  rotational  stiffness  of  the  hypothetical  plate  used  to  calcu¬ 
late  the  reference  frequencies.  The  hypothetical  plate  has  the 
dimensions  of  the  actual  plate  but  with  a  flexural  stiffness 
(El)e(jv  equal  to  the  sum  of  the  flexural  stiffnesses  of  the  upper 
and  lower  face  sheets.  The  rotational  stiffness  at  the  edge  of 
the  hypothetical  plate  is  calculated  for  that  edge  unrestrained 
and  a  clamped  condition  imposed  on  the  opposite  edge. 

Damping  as  a  function  of  the  shear  parameter  for  the  first 
four  modes  of  a  rectangular  sandwich  plate  with  PWU  boundary 
conditions  and  a  geometry  parameter  of  Y  =  3* 5  is  shown  in 
Figures  40  through  42. 

Natural  frequencies  for  sandwich  plates  with  PWU  boundary 
conditions  are  given  in  Figures  43  through  45*  Reference 
frequencies  are  given  in  Table  5* 

A  tabular  presentation  of  the  data  in  Figures  40  through  45 
is  given  in  Tables  24  through  26,  as  well  as  results  for  the 
fifth  and  higher  modes. 
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Figure  39 
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Sandwich  plates  with  PWU  and  PWR 
boundary  conditions 
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Figure  40  Damping  of  a  sandwich  rectangular  plate,  PWU 
boundary  conditions,  Axy  =  1.1,  Y  -  3.5 
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Figure  41  Damping  of  a  sandwich  rectangular  plate,  PWU 
boundary  conditions,  Axy  =2.0,Y=3.5 
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Figure  42  Damping  of  a  sandwich  rectangular  plate ,  PWU 
boundary  conditions,  Axy  =  4.0,  Y  =  3.5 


Figure  43  Natural  frequencies  of  a  rectangular  sandwich  plate 
PWU  boundary  conditions,  Axy  =  1.1,  Y  =  3.5 
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PWU  boundary  conditions ,  Axy 


Natural  frequencies  of  a  rectangular  sandwich  plate 
PWU  boundary  conditions,  Axy  =  4.0,  Y  =  3.5 


TABLE  24 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  pwu  (zero  translation,  elastically  restrained  rotation,  unrestrained  shear) 
Aspect  Ratio  (Axy)  =  1.1 
Geometric  Parameter  (Y)  =  3.5 
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TABLE  25 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PWU  (zero  translation,  elastically  restrained  rotation,  unrestrained  shear) 
Aspect  Ratio  (Axy)  =  2.0 
Geometric  Parameter  (Y)  =  3.5 
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10072.  0.001  10084.  0.001  10340.  0.047  2041.  0.362  2531.  0.308  3118.  0.376 
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3.2.5  PWR  Boundary  Conditions 

Damping  as  a  function  of  the  shear  parameter  for  the  first 
four  modes  of  a  rectangular  sandwich  plate  with  PWR  boundary 
conditions  and  a  geometry  parameter  of  Y  =  3-5  is  shown  in 
Figures  46  through  48. 

Natural  frequencies  for  sandwich  plates  with  PWR  boundary 
conditions  are  shown  in  Figures  49  through  51 •  Reference 
frequencies  are  given  in  Table  5* 

A  tabular  presentation  of  the  data  in  Figures  46  through  51 
is  given  in  Tables  27  through  29,  as  well  as  results  for  the 
fifth  and  higher  modes. 
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Figure  46  Damping  of  a  sandwich  rectangular  plate 
boundary  conditions , -Axy  =  1.1, Y  =  3.5 
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Figure  47  Damping  of  a  sandwich  rectangular  plate,  PWR 
boundary  conditions ,  Axy  =  2.0,  y  =  3*5 
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Figure  48  Damping  of  a  sandwich  rectangular  plate,  PWR 
boundary  conditions,  Axy  =  4.0,  Y  =  3.5 
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Figure  49  Natural  frequencies  of  a  sandwich  rectangular  plate, 
PWR  boundary  conditions,  Axy  =  1.1,  Y  =  3.5 
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Figure  50  Natural  frequencies  of  a  sandwich  rectangular  plate, 
PWR  boundary  conditions,  Axy  =  2.0,  Y  =  3.5 
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Figure  51  Natural  frequencies  of  a  sandwich  rectangular  plate, 
PWR  boundary  conditions,  Axy  =  4.0,  Y  =  3.5 


TABLE  28 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PWR  (zero  translation,  elastically  restrained  rotation,  zero  shear) 
Aspect  Ratio  (Axy)  .  =2.0 

Geometric  Parameter  (Y)  =  3.5 
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10259.  0.001  10324.  0.010  8667.  0.050  2137.  0.416  2591.  0.347  3119.  0.385 


TABLE  29 

MODAL  FREQUENCIES  AND  MODAL  LOSS  FACTORS 

Boundary  Condition  =  PWR  (zero  translation,  elastically  restrained  rotation,  zero  shear) 
Aspect  Ratio  (Axy)  =  4.0 
Geometric  Parameter  (Y)  =  3.5 
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18798.  0.001  18870.  0.007  19286.  0.032  4033.  0.072  4670.  0.238  5744.  0.189 
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3 . 3  EXAMPLE 

As  a  guide  in  using  the  design  charts,  the  sample  problem 
discussed  previously  in  Section  2.4  will  he  solved.  The  follow¬ 
ing  data  are  given: 

Boundary  conditions  =  simply  supported,  unriveted  (PTU) 

Base  layer  thickness,  T-]  =  0.055  inches 

Core  layer  thickness,  T2  =  0.0045  inches 

Constraining  layer  thickness,  T-z  =  0.055  inches 

2 

Viscoelastic  shear  modulus,  Ij2  -  450  Ihf/in 
Base  plate  Young’s  modulus,  E^  =  10  x  10^  Ihf/in^ 
Constraining  layer  Young's  modulus,  E^  =  10  x  10^  Ihf /in^ 
Poisson’s  ratio  of  base  layer, 

=  0.3 

Poisson's  ratio  of  constraining  layer, 

=  0.3 

Mass  density  of  base  layer, 

=0.1  Ibm/in^ 

=  2.59  X  lO""^  Ibf-sec^/in^ 

Mass  density  of  constraining  layer, 

=  0.1  Ibm/in"^ 

=  2.59  X  10“^  Ibf-sec^/in^ 

Mass  density  of  viscoelastic  layer, 

P2  =  0.035  Ibm/in^ 

=  9‘07  X  10"^  Ibf-sec^/in^ 

Core  material  loss  factor,  =  0.3 
Plate  width,  a  =  10  inches 
Plate  length,  b  =  11  inches 

The  dimensionless  variables  that  describe  the  plate  are,  in 
addition  to  n^  already  given: 
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shear  parameter  g  = 


450 


[- 


1 


.0045  -|q7  ^ 

=  40.0 


]11^  X  (1-0.5  ) 


D  =  sum  of  flexural  stiffnesses  of  face  sheets 


5  5 


12(l-vp 

10*^  X  .055^ 
12(1-0.5^) 

=  504.7  Ibf-in 


12(l-v|) 


lo'^  X  .055^ 
12(1-0.5^) 


Y  =  geometry  parameter 
(T^  +  2T2  + 


4D(l-v^) 


[ 


+  E^T^' 


^  (.055  +  2  X  .0045  +  0.55)^  j-io'^  x  .055  x  10*^  x  .05_5_] 
4  X  504.7  X  (1-0.5^)  10*^  X  .055  +  10"^  X  .055 

=  5.51 

Axy  =  in-plane  aspect  ratio 
=  b/a 

=  1 1 .0/10.0 
=  1.10 


The  normalizing  frequency  is  found  from 

p  =  total  plate  mass  per  unit  area 
=  p^T^  +  P2T2  +  ^5^3 

=  2.59xl0“^x.055  +  9.07x10“^x.0045  +  2.59x10“^x.055 
=  2.89x10“^  Ibf-sec^/in^ 
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01 


_  „  1(1)^  +  (1)2-|2 
2Tr  p  ''a''  -• 


1 

2Tr 


1_  D  r  /  If  \  2 

"  P 

304.7 


2,89x10 
=  93.16  Hz 


_  \  (—)^  +  f— )^1 

-5  '■M 1  ^  MO''  ■' 


Figure  6  gives  damping  as  a  function  of  g  for  PTU  boundary 
conditions,  Y  =  3.5,  and  Axy  =  1.1.  Entering  the  chart  with 
g  =  40  gives 


.(1) 


.(2) 


=  0.240 


=  0.336 


for  mode  1 


for  mode  2 


.(5) 


,(4) 


=  0.345 


=  0.331 


for  mode  3 

for  mode  4 


since 

\  =  0.3 

n(') 

=  0.240  X 

0.3  =  0.072 

=  0.336  X 

0.3  =  0.101 

=  0.345  X 

0.3  =  0.104 

=  0.331  X 

0.3  =  0.099 

Natural  frequencies  are  found  by  entering  Figure  21 
(applicable  for  PTU  boundary  conditions,  Y  =  3.5, 
and  Axy  =  1.1)  with  g  =  40  to  obtain 
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1  .75 

5.60 

4.08 

5.63 


Then,  using  the  calculated  reference  frequency  of 
foi  =  93.2  Hz  gives 

f^  =  1 .75  X  93.16  =  1 63.0  Hz 
f2  =  3.60  X  93.16  =  335.4  Hz 
f^  =  4.08  X  93.16  =  580.1  Hz 
f^  =  5.63  X  93.16  =  524.5  Hz 

These  values  may  te  compared  with  results  given  for  the  same 
problem  in  Appendix  A  (the  raw  NASTRAN  output)  and  in  Section  4.3 
(the  closed  form  solution). 
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4.0  CLOSED  FORM  SOLUTION  FOR  HIGH-ORDER  MODES 


4.1  THEORY 

Structures  "built  up  from  plates  always  have  numerous  high- 
order  modes  of  vibration  involving  flexure  of  local  sections. 
Calculation  of  all  these  modes  with  a  single  finite  element 
model,  while  theoretically  possible,  is  neither  practical  nor 
efficient.  However,  such  modes  can  be  important  if  high  fre¬ 
quency  excitation  is  present.  Therefore,  a  method  is  proposed 
for  designing  a  damping  treatment  to  suppress  modes  of  this  type 
which  avoids  the  cost  of  calculating  the  properties  of  a  large 
number  of  essentially  similar  modes.  The  method  is  usable  for 
either  add-on  or  integral  damping  but  is  most  likely  to  be  used 
for  add-on  treatments. 

The  method  is  based  on  the  fact  that  the  higher  order  mode 
shapes  of  rectangular  plates  tend  to  be  sinusoidal  in  both  in¬ 
plane  directions  except  near  the  boundaries.  Boundary  conditions 
have  little  effect  on  the  higher  order  mode  shapes  over  most  of 
the  plate  area.  This  is  true  for  either  classical  uniform  plates 
or  three-layer  sandwich  plates  formed  by  adding  a  constrained 
layer  damping  treatment  to  a  uniform  plate.  This  property  leads 
to  a  useful  relationship  between  natural  frequencies  and  modal 
loss  factors.  When  modal  loss  factors  are  plotted  against  modal 
frequencies  for  a  sandwich  plate,  the  relationship  is  essentially 
independent  of  boundary  conditions  so  long  as  the  boundary 
conditions  themselves  are  non-dissipative  [10]. 

In  the  present  case,  we  are  interested  in  damping  a  number 
of  modes  over  a  fairly  wide  band  of  frequencies.  The  exact  value 
of  natural  frequency  for  each  of  the  many  plate  modes  is  not  of 
particular  importance.  Ve  may  therefore  calculate  the  relation 
between  modal  loss  factor  and  modal  frequency  based  on  any  con¬ 
venient  set  of  boundary  conditions.  The  curve  plotted  for  any 
other  boundary  condition  would  be  formed  by  points  at  different 
frequencies  but  would  still  fall  on  or  near  the  first,  particu¬ 
larly  for  higher  modes.  By  choosing  a  set  of  boundary  conditions 
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which  lead  to  a  simple  closed  form  solution,  we  may  produce  the 
plot  of  damping  vs.  frequency  quite  easily  for  any  given  material 
properties  and  plate  cross-section.  A  few  trials  will  usually  be 
sufficient  to  find  an  appropriate  add-on  treatment  based  on  the 
size,  thickness,  and  material  of  the  base  plate  and  the  desired 
frequency  range. 

The  most  convenient  set  of  boundary  conditions  are  those 
where  all  four  sides  of  the  plate  are  simply  supported  and  shear¬ 
ing  of  the  viscoelastic  core  is  unrestrained.  In  this  case  the 
mode  shapes  are  sinusoidal  all  the  way  to  the  edges  of  the  plate. 
Two  closed  form  solutions  for  this  case  are  available  in  the 
literature  [8,10].  The  former  solution,  due  to  Abdulhadi,  has 
been  used  in  this  work  because  the  latter  could  not  be  made  to 
produce  results  in  agreement  with  MSC/NASTRAN-MSE  even  for  small 
values  of  the  core  material  loss  factor.  The  Abdulhadi  solution 
did  produce  good  agreement  for  core  loss  factors  small  compared 
to  unity.  It  diverged  somewhat  as  the  core  loss  factor 
approached  unity,  as  shown  in  Section  2  of  this  report  and  in 
previous  work  by  the  authors  [1]. 

The  closed  form  solution  from  Ref.  [8]  for  the  natural  fre¬ 
quency  and  modal  loss  factor  of  a  simply  supported  rectangular 
sandwich  plate  is,  using  the  notation  of  Eq.  (6-8): 


6^K  T^K  ct^  P 

e  f  2  e  r  ,  ^ 

2  -ir—  h - +  1  +  »v) 

a  D  (jo 

=  TpKga^  2  2  ^ 

^2 


(28) 


(29) 
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where 


a  2  _ 

=  (-^  + 1 )  +4 


(50) 


^  (1  +  +  4^) 

®2 


(51 ) 


i;,t5 


12(1-vp 


"3^3 

12(1-v2) 


(52) 


T  T 
1  *=5 


(53) 


p  =  ^1*^1  ^2^2  ^3^3 


(54) 


(1+v,  )  [E,T,  (l-v_)  +  E-T-(l-vJ] 


(55) 


=  radian  frequency  of  the  r*th  mode 
=  structural  loss  factor  of  the  r'th  mode 
=  material  loss  factor  of  the  viscoelastic  core 
=  thickness  of  the  core  layer 


Ti,T3 


''l 


=  thicknesses  of  the  face  sheets 


E^  » ^2  >  ~ 


Young’s  moduli  of  the  materials  of  the 
three  layers 


=  Poisson’s  ratios  for  the  three  layers 


=  complex  shear  modulus  of  the  core 

=  ifgd  +  i  r,^) 


,mTt»2  /nir\2 

(^)  +  (^) 
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ajb  =  in-plane  dimensions  of  the  plate 

m,n  =  integers 

The  use  of  simply  supported  plate  solutions  for  designing  a 
damping  treatment  for  a  plate  with  any  non-dissipative  boundary 
conditions  may  seem  like  a  gross  approximation.  However,  the  end 
purpose  of  the  damping  must  be  kept  in  mind.  Local  plate  modes 
of  an  all-welded  structure  without  any  damping  treatment  would 
typically  have  loss  factors  on  the  order  of  .001.  Loss  factors 
predicted  by  the  simple  method  described  here  are  on  the  order  of 
.1,  even  without  any  extensive  searching  for  an  optimum  damping 
treatment.  If  boundary  condition  effects  change  the  loss  factors 
by  a  factor  of  three,  the  end  conclusion  regarding  the  damping 
treatment  will  remain  the  same;  namely,  that  it  produces  a  sub¬ 
stantial  (factor  of  30)  reduction  in  resonant  response  to 
periodic  input  forces.  Complex  eigenvalue  solutions  for  sandwich 
beams  have  shown  that  the  variation  of  damping  with  boundary 
conditions  is  small  for  higher  modes  and  generally  less  than  a 
factor  of  three  for  low  order  modes  [l]. 

4.2  SOFTWARE  IMPLEMENTATION 

The  eq^uations  given  above  have  been  implemented  in  an  inter¬ 
active  program  called  SPLT61  [11  ]•  The  program  allows  a  designer 
to  quickly  evaluate  a  number  of  possible  constrained  layer 
treatments  with  negligible  cost  for  computing.  The  input  to  the 
program  includes  flexural  wavelength  in  each  direction.  These 
are  usually  set  equal  to  twice  the  plate  dimensions,  i.e.,  the 
exact  value  for  a  simply  supported  plate.  The  output  is  the 
exact  (i.e.,  closed  form,  complex  eigenvalue)  solution  for  modal 
frequencies  and  loss  factors  of  a  simply  supported  sandwich 
plate.  As  noted  above,  this  frequency-damping  relationship  is 
also  correct  in  the  limit  of  increasing  mode  number  for  other 
boundary  conditions. 
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The  validity  of  using  a  solution  for  a  simply  supported 
sandwich  plate  to  predict  damping  for  other  boundary  conditions 
was  tested  as  follows.  The  SPLT61  program  was  used  to  obtain 
modal  damping  as  a  function  of  modal  frequency  for  a  number  of 
different  values  of  the  shear  parameter.  The  shear  parameter  was 
varied  by  changing  G2,  the  core  shear  modulus.  The  SPLT61 
results  represented  an  exact  solution  for  a  simply  supported 
plate,  (i.e.,  PTU  boundary  conditions).  A  similar  analysis  was 
performed  for  a  plate  with  PLR  (fixed)  boundary  conditions  using 
the  modal  strain  energy  method.  A  comparison  of  results  is  shown 
in  Figures  52  through  55*  The  damping  parameter  is  plotted  vs. 
normalized  modal  frequency  for  various  values  of  the  shear  para¬ 
meter.  The  frequency  is  normalized  on  a  reference  frequency  fQ^ 
as  described  in  a  previous  section.  The  same  information  is 
shown  in  dimensional  form  in  Figures  56  through  60.  It  may  be 
seen  that  the  results  for  the  drastically  different  boundary  con¬ 
ditions  do,  in  fact,  converge  at  high  frequency.  Significantly, 
the  rate  of  convergence  depends  on  the  shear  parameter. 

4.5  EXAMPLE 

The  use  of  the  SPLT61  program  is  illustrated  below.  The 
physical  situation  being  analyzed  is  the  same  one  treated  using 
the  NASTRAN  modal  strain  energy  method  in  Section  2.4  and  using 
the  design  charts  in  Section  5*3*  Since  true  simply  supported 
boundary  conditions  are  assumed,  and  the  viscoelastic  loss  factor 
is  small,  the  SPLT61  solution  agrees  closely  with  the  NASTRAN/MSE 
results . 

The  program  begins  by  displaying  default  values  for  all  the 
physical  parameters  describing  a  three-layer,  rectangular 
sandwich  plate.  The  user  is  prompted  to  change  any  or  all  of  the 
values,  to  store  the  entire  list  in  a  disc  file,  or  to  proceed 
with  the  calculations.  He  has  the  option  of  resetting  the  entire 
list  of  values  to  a  set  previously  stored  on  disc.  Once  the  user 
is  satisfied  with  the  input  data,  he  commands  the  program  to 
proceed.  It  then  calculates  a  table  of  natural  frequencies  and 
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modal  loss  factors  for  the  first  ten  allowahle  values  of  wave- 
number  in  each  in-plane  direction  (i.e.,  one  hundred  normal 
modes).  The  resulting  table,  along  with  the  input  parameter 
list,  is  written  to  a  disc  file  in  ASCII  format  to  be  printed 
later.  The  modes  calculated  are  not  necessarily  the  lowest  one 
hundred,  but  will  always  contain  the  lowest  ten. 

The  output  file  for  the  sample  case  is  shown  in  Table  30. 

The  loss  factors  are  given  directly  rather  than  being  normalized 
on  the  core  material  loss  factor.  They  must  therefore  be  divided 
by  that  value  (0.3  for  the  sample  case)  for  comparison  with 
results  from  the  other  two  methods. 
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Effect  of  boundary  conditions  on  damping  for 
rectangular  sandwich  plates,  g  =  1,3  and  10 
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Figure  53  Effect  of  boundary  conditions  on  damping  for 
rectangular  sandwich  plates,  g  =  40 
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Figure  56  Effect  of  boundary  conditions  on  damping  for 
rectangular  sandwich  plates,  g  =  1,3  and  10 
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Figure  57  Effect  of  boundary  conditions  on  damping  for 
rectangular  sandwich  plates,  g  =  40 
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Effect  of  boundary  conditions  on  damping  for 
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Aerospace  Structures 
Information  and  Analysis  Center 


c 

D- 

OH 

c 

-J 

rD 

o 

c 

I— 

o 

UJ 

q: 

o 

LjJ 

l-“ 

DC. 

O 

a. 

D_ 

rD 

CO 

>- 

D. 


O 


O 

CO 


q: 

o 


CO 

o 


o 

DC 


CO 


13 

CO 


o 

CO 


CO 

c 


o 

►- 

o 


Q. 

Ol 


ru 

rsx 


Ol 

fn  nj 


•  CO 

tc  to 
tn  ru 
in  • 
in 


111 

•  ir 

^  nj 
• 
tn 


m  4> 
^  ru 


o  r\j  in  oj 

CT'  :•  ^  • 


ft 

01 

41 

•  m 

o 

9-  m 


t  •  f  • 

lU  UJ  Uj  Ui  Ul  UJ 

c;9OC><90VO0000 

OOOOOOOO  OOOOOO 

o  O  O  O  '*0  oooooooo 
mminooooooinoo  •  • 
in^in^tn»^rnrn»^»ni>«moru- 

»  •  •  •  a-  ••••••  •oj  nj 

ooo=>  ooooooo 


*-»  9> 

X  to 
1-  UJ  to  >» 
^  *o 

a:  UJ  X 

MJ  w  ^  Ui 

>>  M  :ic  >• 

X  o  < 


*  * 

• 

•  » 

•  • 

• 

9 

•  • 

•  • 

9 

• 

•  • 

•  • 

9 

• 

•  • 

»  • 

m 

•  » 

•  m 

• 

•  » 

m  9 

9 

•  • 

9  9 

9 

• 

•  • 

9  9 

9 

» 

•  • 

9  9 

9 

• 

•  • 

9 '  m 

•  X 

9  9 

• 

•  UJ 

9  9 

m 

• 

•  >- 

•  X  X 

9 

• 

•  c 

•  UJ  UJ 

9 

o 

• 

•  J 

• 

>-  >- 

9 

U|  UJ 

-• 


3 

3 


o  O  ►-  t 
z  *-<  z  *■ 
*-H  •-  CC  •-»  * 
z  V7  'Ui  iCt 
l-l  4*  >-  w  • 

<  -J  «  <  - 
X  ;aJ  J  X 
1-0  ►- 
to  U<UI  to 
Z  to  CO  z 

3  *-i  <  o 

o  >  ao  o 


>  to  to 

ID  O 

>  ^  -I 

:  3 

)  o  o 
o  o 
X  X 

JX  to 

4  • 

Ui  C9 

r  z 

to  D 

O 


O  •  .1  .1 

z  • 

1-4  X  C9  O 
z  UI  z  ^ 

»-i  >-  •-<  i—  X 

<r  <  z  to  UJ 

X  -j  •-•  «*  >- 
*-  4X  3 
CO  UJ  X  UJ 
Z  CO  •—  O 

o  <t  to  o 

O  Si  Z  CO 

a  -• 

u.  U.  U  » 

o  o 


u  t_> 
UJ  UJ 
X  X 


•  o  o 

_l  X 

O  X  >“ 

ui  t~ 

<0  U>  tZ  z 

«c  <  ^ 

cu  u. 

X  X 


•  t^ 

in  ^ 
in  m 


»  «4  •  (D 


m  sr 
so  to 
in  • 


r  nj 

t^  m 


tr 

o  • 

tn 


•  m 

-4  O' 
r*-  f\i 
St  • 

in 


ru 


I 

UJ 

in 


•  r\l  •  *4 


4> 

m  rvi 
a*  ■• 

»n  : 


a  •  44 


CD  -4 

fo  rvj 
so  • 


UJ 

m 

•  in 


tn 

•  tn 

o 


ft 

Ui 


•  o* 

m  4* 

aD  9 


-4  » 
t—  • 

m 


•  o* 

fO  cp 
nj  to 


•  o 

-4  a 

O  fO 


•  in 

«  ID 
-4  r\j 

4>  • 

in 


•  o 

44  a 

U>  *M 
f\l  • 


199 

91^ 

o 

o 

o 

o 

o 

o 

■  o 

O 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

til 

Ui 

Ul 

UJ 

UJ 

Ui 

UJ 

Ul 

0 

9 

m 

0 

9 

ru 

•  ao 

•  m 

•  ru 

•  r«- 

•  o 

•  in 

•  ru 

44  o 

a  a 

ru  9 

a  tn 

fu  a 

0  9 

0  ru 

m  a 

44  in 

•n  9 

in  9 

9  9 

m  m 

a  ro 

9  m 

0  ru 

0  • 

a  • 

44  • 

9  • 

0  • 

ru  :• 

0  • 

9  • 

ru 

ru 

»o 

tn  ■ 

9  ; 

9 

in 

•  in 
a  tn 
4>  rv 
O'  • 
D 


ft 

UJ 


•4  CO 


ft 

uJ 

<c 


in  m 
o  ru 


ft 

Ui 


4> 
•  4> 
O- 
^  in 
ru  • 
ru 


1  I 


ft 

Ui 


44  in 

eo  U^ 


•  «■ 
m  -D 
9  r 

44  4>  •  44 

to  '  to 


I 

UJ 
o 
•  O' 

ru  4> 

CO  tn 


ft 

UJ 
D 
•  O 

rw  0 

m  ru 
in  • 
in 


•  iO 

o  m 

a  Ai 

m  • 


to  « 

in  '• 


ft 

UJ 

4> 

•  o 
D-  ru 
a  r»». 

D  • 

44  *4  ru 


i 

Ui 


o  r«- 

.4  0 

o 


<X>  ST 

9  in 

in  • 


•  *n 

lA  ^ 
0  a 


I 

ul 
ID 
•  m 
o  9 
m  tn 
ru  • 


Ol 

ru 

•  m 

tn  44 

ao  O' 

o  • 


o 

o 

o 

O 

o 

o 

o 

• 

« 

• 

• 

t 

• 

• 

Oj 

iZl 

Ol 

iXJ 

Ul 

lU 

i£j 

0 

44 

9 

r* 

o 

0 

•  0 

•  r< 

•  ru 

•  m 

•  a 

•  in 

•  o 

f'-  9 

14.  a 

ru  -4 

o  ru 

44  Kl 

tn  9 

ru  o 

•4  0 

ro 

9 

tn  9 

o  in 

in  9 

a  9 

ru  • 

9  • 

r*  • 

44  • 

9  • 

44  • 

r*.  • 

44 

ru  - 

ru 

tn 

m 

44 

44 

44 

O 

o 

o 

o 

O 

a 

1 

ft 

ft 

ft 

ft 

ft 

UJ 

UJ 

Ul 

UJ 

UJ 

ul 

*4 

9 

n* 

tn 

44  ■ 

44 

0 

0  O 

•  ru 

• 

•  9 

•  ru 

•  44 

•'9 

•  o 

ru  ru 

tn  -* 

m  a 

0  o 

9  44 

fU  rO 

m  44 

9  a 

rw  0 

9 

f-l  0 

a  m 

fo  a 

m  • 

o  • 

in  ♦ 

1^  *• 

ru  • 

P4  • 

9  *• 

0 

•* 

ru 

ru 

fO 

o 

o 

o 

o 

o 

1 

1 

1 

ft 

ft 

Ul 

Ul 

Ul 

UJ 

UJ 

•  a> 

a  44 

r*.  m 


tli 

9  ■ 

•  m 


fo  a  in  ^ 


to  44 

m  tn 

9  '• 


O'  44 

•  o 

ST  — 
CO  • 


•  r^. 

9  O 
O  9 

44  •  tn 


•  to 
ru  rw 


•  9 

9  sr 

»0  u» 

m  • 


in  9 

9 


r*  O' 

9  tn 


44  •  0 


44  -4  ,  44  ru 


i  s 


u.  u.  u.  to  • 


44 


in  ’•  44  • 


I 

UJ< 

ru 


•  -o 
o  a 


9  **  Kl  • 


I  ft 

ul  UJ> 

ru  44 

•  ?••  •  ru 

1  ru  ru  a 

lU  9 


•  9 

44 

«  U^ 


44  44  0  »0  O 


a 

9 

in 


-4  :  44  ru 


o  O  o  o 


3 

O  UJ 

to  to 

>4*  < 


X  X  44 

to  to 
I  0  CO  z  z 
>  C  S  UJ  UJ 

:  z  z  o  o 

t  3  O 

.  44  ►-*  I—  »— 

0  0  X  X 

0  0  3  3 


O  O  UJ  UJ 
X  X  « 


0  0  3  3 
O  Z  Z 
>>  3  UJ  Ul 
►4  -I  -I 

44  (_>  UJ  UJ  I 

0  44  >  >• 

Z  ♦-  <  « 
UJ  0  £  Z 
Q  4t 

^  X  X 
t-  UJ  3  3 
r  o  X  X 
3  Ui'-t  *-« 

44  0  X  X 
Jj  44  4C  < 

«  <  >  X  X 


X 

3  ‘X 
z  •-* 
uJ  O 
-i 

UllX 


44-  ^  tn 


44  ru 


9 

9  -4 


m  *0  ru 


-O 


121 


3a3<»t*0l  .333<>E*0l 


Aerospace  Structures 
Information  and  Analysis  Center 


5.0  SUMMARY  AND  CONCLUSIONS 

Three  methods  have  been  presented  for  the  dynamic  analysis 
and  design  of  viscoelastically  damped  sandwich  plates.  The 
methods  are  complementary  in  that  each  represents  a  different 
trade-off  of  accuracy,  generality,  and  cost  of  use.  The 
theoretical  basis  of  each  has  been  described  along  with  sample 
problems . 

The  most  general  of  the  three  is  the  modal  strain  energy 
method  implemented  in  MSC/NASTRAN.  It  is  fairly  new,  having  been 
developed  primarily  by  the  authors  of  this  report.  It  can  accom¬ 
modate  virtually  any  combination  of  plate  geometry  and  boundary 
conditions.  The  method  is  by  no  means  limited  to  sandwich 
plates,  although  it  has  seen  extensive  application  there  due  to 
the  efficiency  of  this  construction  for  vibration  damping.  The 
price  for  the  accuracy  and  generality  of  MSE  is  that  the  user 
must  be  fluent  in  NASTRAN  and  must  prepare  and  run  a  finite 
element  model  for  each  candidate  design. 

A  simplified  version  of  MSC/NASTRAN-MSE  for  rectangular 
sandwich  plates  is  presented  in  the  form  of  design  charts  derived 
from  a  large  number  of  MSC/NASTRAN-MSE  runs.  They  allow  the  user 
to  rapidly  obtain  values  for  modal  loss  factors  and  modal  fre¬ 
quencies  of  sandwich  plates  with  various  boundary  conditions 
Plate  geometry  and  material  properties  are  specified  in  terms  of 
dimensionless  groups  to  allow  the  maximum  information  to  be 
conveyed  by  each  chart. 

The  simplest  of  the  three  methods  is  based  on  the  use  of  a 
closed  form  solution  that  is  strictly  applicable  only  to  simply 
supported  rectangular  sandwich  plates.  It  is  shown  that  the 
solution  may  be  used  with  other  boundary  conditions  to  obtain 
damping  estimates  of  useful  accuracy  for  higher  order  modes.  The 
method  involves  negligible  costs  for  computation  and  has  been 
implemented  in  an  interactive  Fortran  program. 

The  latter  two  methods  are  applicable  only  to  single  rectan¬ 
gular  plates  rather  than  to  assemblages  built  up  of  plate 
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elements.  They  are  nonetheless  useful  in  that  a  designer  often 
seeks  to  increase  the  damping  of  local  modes  of  individual  plate 
sections.  The  modal  strain  energy  method,  when  implemented  in 
NASTRAN,  is  quite  general  and  will  readily  accommodate  huilt-up 
structures  with  integral  damping. 
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APPENDIX  A 

SAMPLE  INPUT  AND  OUTPUT  FOR 
NASTRAN/MODAL  STRAIN  ENERGY  ANALYSIS 
OF  A  SANDWICH  PLATE 
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Sample  NASTRAN  input  is  given  for  the  following  case: 

Boundary  condition  =  PTU  (zero  out-of-plane  translation, 

zero  moment,  unrestrained  core  shear) 

Pace  sheet  thicknesses  (equal),  =  0.055  in. 

Core  layer  thickness,  T2  =  0.0045  in. 

2 

Viscoelastic  shear  modulus,  Ir^  =  450  Ibf/in 

7  2 

Face  sheet  Young's  moduli  (equal),  =  10'  Ibf/in 

Poisson's  ratio  of  face  sheets  (equal),  =  0.3 

Poisson's  ratio  of  core  layer,  =  0.49990 
Mass  density  of  face  sheets  (equal), 

=  2.59  X  10“^  Ibf-sec^/in^ 
2/4 

Mass  density  of  core  layer,  p^  =  9.07  x  10  Ibf-sec  /in 
Plate  in-plane  dimensions,  a  x  b  =  10  x  11  inches 


The  entire  bulk  data  deck  is  listed  for  clarity  although 
most  of  it  was  produced  automatically  by  a  mesh  generator. 
Executive  and  case  control  decks  are  also  listed.  The  grid  and 
the  numbering  system  is  illustrated  in  Figure  A-1 . 

Output  listed  for  the  sample  case  includes  the  first  page  of 
the  eigenvalue  table,  mode  shapes  for  the  first  four  modes,  and 
strain  energy  distributions  for  the  first  four  modes.  Set  99 
includes  all  the  solid  elements  used  to  model  the  core.  Thus, 
the  "percent  of  total"  figure  printed  out  is  (after  dividing  by 
100)  exactly  the  strain  energy  fraction  of  Eq.  4,  which 

equates  to  the  loss  parameter 


A- 2 


1  101  201  301  401  501  601  701  801  901  1001  1101 


101  201  301  401  501  601  701  801  901  1001  1101 


Figure  A-1 


Finite  element  grid  for  modal  striain  energy 
analysis  of  rectangular  sandwich^plate.  Upper 
face  sheet  showing  partial  grid  numbering  and 
QUAD4  element  numbering. 
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Figure  A- 2  Finite  element  grid  for  modal  strain  energy 

analysis  of  rectangular  sandwich  plate. 
Viscoelastic  core  layer  showing  partial 
HEXA  element  numbering. 
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Figure  A- 3 


Finite  element  grid  for  modal  strain  energy 
analysis  of  rectangular  sandwich  plate. 
Lower  face  sheet  showing  grid  numbering 
and  QUAD4  element  numbering. 
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